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Abstract 

In conventional approaches to the homogenization of random particulate composites, both 
the distribution and size of the component phase particles are often inadequately taken into 
account. Commonly, the spatial distributions are characterized by volume fraction alone, 
while the electromagnetic response of each component particle is represented as a vanishingly 
small depolarization volume. The strong-permittivity-fluctuation theory (SPFT) provides 
an alternative approach to homogenization wherein a comprehensive description of distri- 
butional statistics of the component phases is accommodated. The bilocally-approximated 
SPFT is presented here for the anisotropic homogenized composite which arises from com- 
ponent phases comprising ellipsoidal particles. The distribution of the component phases is 
characterized by a two-point correlation function and its associated correlation length. Each 
component phase particle is represented as an ellipsoidal depolarization region of nonzero 
volume. The effects of depolarization volume and correlation length are investigated through 
considering representative numerical examples. It is demonstrated that both the spatial ex- 
tent of the component phase particles and their spatial distributions are important factors 
in estimating coherent scattering losses of the macroscopic field. 
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1 Introduction 



Consider the propagation of electromagnetic radiation through a composite medium com- 
prising a random distribution of particles. Provided wavelengths are sufficiently long com- 
pared with the dimensions of the component particles, the composite may be regarded as 
an effectively homogenous medium. The estimation of the constitutive parameters of such a 
homogenized composite medium (HCM) is a matter of long-standing scientific and technolog- 
ical importance [1. 2]. Furthermore, recent advances relating to HCM-based metamaterials 
have promoted interest in this area and highlighted the necessity for accurate formalisms to 
estimate the constitutive parameters of complex HCMs [3]. 

The microstructural details of the component phases are often inadequately incorporated 
in homogenization formalisms. In particular, conventional approaches to homogenization, 
as exemplified by the widely-applied Maxwell Garnett and Bruggeman formalisms, gener- 
ally involve simplistic treatments of the distributional statistics and sizes of the compo- 
nent phase particles [4]. It is noteworthy that both the Maxwell-Garnett and Bruggeman 
formalisms arise in the long-wavelength regime. The long-wavelength derivation of the 
Maxwell-Garnett formalism follows a from rigorous treatment of the singularity of the free- 
space dyadic Green function [5, 6]. The Bruggeman formalism emerges from the strong- 
permittivity- fluctuation theory (SPFT) under the long-wavelength approximation [7, 8]. 

The SPFT provides an alternative approach to homogenization in which a comprehensive de- 
scription of the distributional statistics of the component phases is accommodated. Though 
the SPFT was originally developed for wave propagation in continuous random mediums 
[7, 8, 9, 10], it has more recently gained prominence in the homogenization of particulate 
composites [11, 12, 13, 14, 15]. Within the SPFT, estimates of the HCM constitutive param- 
eters are calculated as successive refinements to the constitutive parameters of a homogenous 
comparison medium. Iterates are expressed in terms of correlation functions describing the 
spatial distribution of the component phases. Correlation functions of arbitrarily high order 
may be incorporated in principle, but the SPFT is usually implemented at the second order 
level, known as the bilocal approximation. A two-point correlation function and its associ- 
ated correlation length characterize the component phase distributions for the second order 
SPFT. At lowest order (i.e., zeroth and first order), the SPFT estimate of HCM permittivity 
is identical to that of the Bruggeman homogenization formalism [2] . 

Depolarization dyadics are key ingredients in homogenization analyses as they character- 
ize the electromagnetic field inside inclusions (i.e., component phase particles) embedded 
within a homogenous background. Commonly, the inclusion particles are treated as vanish- 
ingly small point-like entities; under this approximation, the corresponding depolarization 
dyadic is represented by the singularity of the associated dyadic Green function [16, 17]. 
Through neglecting the spatial extent of the inclusions, potentially important information 
is lost, particularly when coherent scattering losses are to be considered [18, 19]. It is noted 
that extended versions of both the Maxwell Garnett formalism [20, 21] and the Bruggeman 
formalism [20, 22] have been developed in which a nonzero volume is attributed to the com- 
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ponent phase particles. However, these analyses apply only to isotropic HCMs and adopt a 
simplistic description of the distributional statistics of the component phases. 

In the present study, the SPFT is presented for an anisotropic dielectric HCM. The analysis 
takes place within the frequency-domain, wherein dissipation is signalled by constitutive 
parameters having nonzero imaginary parts. The component phases are taken to consist of 
randomly distributed, ellipsoidal particles of nonzero volume. In Section 2, the depolarization 
dyadic appropriate to an ellipsoidal inclusion of linear dimensions proportional to the size 
parameter 77 > 0, embedded within an anisotropic comparison medium, is developed. Details 
of the SPFT-based homogenization are presented in Section 3. In Section 4, the theoretical 
results are illustrated by representative numerical examples relating to both nondissipative 
and dissipative component phases. In particular, nondissipative component phases giving 
rise to dissipative HCMs — through radiative scattering loss from the macroscopic coherent 
field [23] — is highhghted. A few concluding remarks are provided in Section 5. 

In the notation adopted, vector quantities are underlined whereas 3x3 dyadics are double 
underlined. The inverse, adjoint, determinant and trace of a dyadic M are denoted by M 
M"* , det [M] andtr [ M ], respectively. The identity dyadic is represented by /. All field- 
related quantities are implicitly functions of the angular frquency uo. The permittivity and 
permeability of free space arc denoted as e,, and yUo, respectively; the free-space wavenumbcr 
is ko = ujy/e^. The real and imaginary parts of complex- valued quantities are represented 
by Re ( • ) and Im ( • ) , respectively. 



2 Depolarization 

Let us consider an ellipsoidal inclusion of volume VJ', orientated with its principal axes lying 
along the Cartesian x, y, z axes. The ellipsoidal surface of is parameterized by 

uoA)^vg'm<i>), (1) 

where f {9, 0) is the radial unit vector specified by the spherical polar coordinates 9 and 4>. 
The diagonal shape dyadic 

[/ = -^=4== diag(C/„ {U,,Uy,U,>0), (2) 

maps the spherical region of radius 77 onto the ellipsoidal region VJ'. The linear dimensions 
of the ellipsoidal inclusion, as determined by 77, are assumed to be sufficiently small that the 
electromagnetic long-wavelength regime pertains, but not vanishingly small. 

Suppose the ellipsoidal inclusion is embedded within a homogenous comparison medium. In 
consonance with the shape dyadic (2), the comparison medium is taken as an anisotropic 
dielectric medium characterized by the diagonal permittivity dyadic e with principal axes 
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aligned with those of IJ_. We remark that e^^, being estabhshed in the frequency-domain, is 
imphcitly a fuction of u. The electromagnetic response of the ellipsoidal inclusion is provided 
by the depolarization dyadic 

D = [ G ir)d^r=[ G {U'r)d^r. (3) 

Herein, G^{r) is the dyadic Green function of the comparison medium which satisfies the 
nonhomogenous vector Helmholtz equation 

(vxVxI-^VoiJ 'g^Jr-r:) = iuJiioHr-i::)L (4) 

Explicit representations of Green functions are not generally available for anisotropic medi- 
ums [24]. However, it suffices for our present purposes to consider the Fourier transform of 
G^{r), namely 

£cm(^) ^ l=cm^-^ exp(-ig • r) d^r. (5) 
Taking the Fourier transform of equation (4) delivers G (q) as 

Thereby, equation (3) yields [16, 17] 

In order to consider the depolarization of an inclusion of nonzero volume, we express D as 
the sum 

D^DJ' + D^ (8) 



where 



with 







= 27r2 , 








= 27r2 . 




= cm 






•1) 



qr] 



- cosiqf])^ ^'^ig-' ' q) d\, (9) 
-cos{qrj)^ a'^Jg-'-q)d% (10) 



aj£-^-g)-gI(r^-t)' (11) 



lim G -q). (12) 



g— >oo 



Thus, the dyadic represents the depolarization contribution arising from the vanishingly 
small region limVT', whereas the dyadic D'' provides the depolarization contribution arising 



from the region of nonzero volume ( — limVJ' 1 . In homogenization studies, it is common 

practice to neglect _D'' and assume that the depolarization dyadic is given by Z2° alone [2]. 
However, studies of isotropic HCMs have emphasized the importance of the nonzero spatial 
extent of depolarization regions [18, 19, 20, 21, 22]. 

The properties of depolarization dyadics associated with vanishingly small inclusions have 
been widely investigated: the volume integral (10) simplifies to the 77-independent surface 
integral [16, 17] 

= / dcj) I d0 sm0 I , ^ I , (13) 



V tr e 



(e -a) 

\ =cm = / 



where 



/ sin^ 9 cos^ 6 sin^ 9 sin^ 6 cos^ ^ \ , , 

4^d,ag(^^.^^.^). (14) 

Furthermore, the integrations of (13) reduce to elliptic function representations [25]. In the 
case of spheroidal inclusion geometries, hyperbolic functions provide an evaluation of [16], 
while for the degenerate isotropic case — Uy — the well-known result ^ — (l/3ia;) 
emerges [26]. 

We focus attention on . By analogy with an equivalent integration which arises within 
the strong-permittivity-fluctuation theory [27], the application of residue calculus to (9) 
provides 

£' = ^m^)^ (15) 

where the dyadic function W_{r]) has the surface integral representation 

(16) 
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with 



a 



2 e - tr ( e ) / 

cm V =cm 



, , , , tr f •^') - ftr f e'^'^'O tr (A) 

\ = cm = / 



/3 = e' 



A 



det e 

ad; \=c"^ 

tr ( e -A 

cm 



A: 



trf e"'^^ 'A)-ti( e"*) tr(A) 

\ = cm ^ / \ = cm / V ^ / 



- 4 det ( e ) tr ( A ) tr f e -A 

_ = cm / v = / \^ = cm = 



tr f e""^-?' ' 

\ = cm ) 


tr(4) 


- tr f e"''^' 

y =cm 


•4) 









(17) 
(18) 

(19) 
(20) 



The surface integrals of (13) and (16) are straightforwardly evaluated by standard numerical 
techniques [28]. 



3 Homogenization 



The SPFT may be applied to estimate the constitutive parameters of HCMs. Let us consider 
the homogenization of a two-phase composite wherein the two component phases, labelled as 
a and 6, comprise ellipsoidal particles of shape specified by C/ and linear dimensions specified 
by ?7 > 0. A random distribution of identically orientated particles is envisaged. The 
component phases a and h are taken to be isotropic dielectric mediums with permittivities 
and 65, respectively. 

The distributional statistics of the component phases are described in terms of moments of 
the characteristic functions 



$^(r) 



1, 
0, 



{i = a,b). 



(21) 



The volume fraction of phase £, namely fg , is given by the first statistical moment of $^ ; 
i.e., {^i{r) ) = fe ■ Clearly, fa + fb = 1- The second statistical moment of $^ provides a 
two-point covariance function. We adopt the physically-motivated form [29] 



{MD Mr:). 



Mr.)) {Mr.')), 
Mr)): 



Ur^ ' ir 
c/-^ • (r 



r!)\> L, 
r!)\<L, 



(22) 
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where L > is the correlation length. It is remarked that the specific form of the covariance 
function has httle influence on SPFT estimates of HCM constititutive parameters [30] . 

The nth order SPFT estimate of the HCM permittivity, namely e'"^, is based upon the 
iterative refinement of the comparison medium permittivity, namely e . To zeroth order 
and first order, the SPFT permittivity estimate is identical to the comparison medium 
permittivity [2]; i.e., 

= elij = e . (23) 

= ab =ao =cm ^ ' 

The well-known Bruggeman homogenization formalism provides the estimate of e [14]. 
Thereby, e emerges through solving the nonlinear equations 



/aX +/6X =0, 

= a =0 ~ 



wherein the polarization dyadics 



X 



(£ = a,6). 



are dependent upon the inclusion size parameter r\ via 



I \ w D ' { til - t \ 

= = \ = =cm / 

Recursive procedures generate the pth iterate of e as 

° =cm 

lc>] = ^{lc.b-i]}, 

with the operator T being defined by 



(24) 



(25) 



(26) 



(27) 



(28) 



The SPFT is most widely implemented at the second order level — known as the bilocal 
approximation — which provides the following estimate of the HCM permittivity dyadic [14] 



= a6 =( 



iu ^- — — ' 



E. 



(29) 



Thus, the inclusion size parameter 77 infiuences e'^j^ directly through the depolarization dyadic 



and indirectly through the mass operator [9] dyadic term 



fafb 



(30) 



Notice that the correlation length L — which plays a key role in the second order SPFT — is 
not relevant to the zeroth order SPFT. 
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4 Numerical studies 



Wc investigate the theoretical results presented in Sections 2 and 3 by means of represen- 
tative numerical studies. To highlight the particular effects of depolarization volume and 
correlation length, it is expedient to begin with a study of HCMs arising from nondissipa- 
tive component mediums (i.e., £„,(, G M) before proceeding to consider HCMs arising from 
dissipative component mediums (i.e., €afi £ C with Im (ea^) > 0). However, we note that 
caution should be exercised when applying the SPFT and Bruggeman formalism to weakly 
dissipative HCMs where Re (e^) and Re (e;,) have opposite signs [31]. For all calculations 
reported here, the angular frequency u) is fixed at 2% x 10^° rad s~^. 

4.1 Nondissipative component phases 

Consider, for example, the HCM which arises from the nondissipative component phases 
specified by the permittivities and shape parameters 

€a = 2eo, 6fe = 12eo 

= 1, Uy = 3, U, = 15 

In Figure 1, the components of the zeroth order SPFT estimate of HCM permittivity 
e'^l^ = eodiag(e^^Q, e^^Q, e^^g) are plotted against volume fraction and relative inclusion 
size koTj. The permittivity parameters e^f^^ are constrained to agree with ej, and at /a = 
and fa = 1, respectively. For volume fractions G (0, 1), the ellipsoidal particulate ge- 
ometry of the component phases endows the HCM with biaxial anisotropy and we observe 
that e^(,Q, e^^Q, e„^Q, and each take different values. At ?7 = the estimates of the HCM 
permittivity parameters ej^^g'^ are the same as those of the conventional Bruggeman homoge- 
nization formalism. When rj > the parameters e^f^^ become complex-valued with positive 
imaginary parts, despite the component phase permittivities and being real-vahied. 
The manifestation of Im (e^^g'^) > indicates radiative scattering loss from the macroscopic 
coherent field [23]. Thus, we see that by attributing a nonzero volume to the component 
phase inclusions, coherent scattering is accommodated, thereby giving rise to a dissipative 
HCM. Since the imaginary parts of e^^^'^ are observed to increase with increasing rj, it is 
deduced that more radiation is scattered out of the coherent field by enlarging the depolar- 
ization volumes. The real parts of the HCM permittivity parameters e^w^ found to be 
relatively insensitive to the inclusion size parameter rj. 

Let us turn now to the bilocally-approximated SPFT. For vanishingly small component 
phase inclusions (i.e., 77 = 0), the components of the second order SPFT estimate of HCM 
permittivity e^^j^ = €0 diag (e^^g) ^ab2' ^062) graphed against volume fraction and relative 
correlation length koL in Figure 2. As L — > 0, the second order SPFT estimates converge 
to those of the zeroth order SPFT. Accordingly, the estimates of the HCM permittivity 
parameters e^^2^ at L = are the same as the estimates deriving from the conventional 
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Bruggeman homogcnization formalism. As the correlation length L increases from 0, the 

HCM parameters ^^^^2^ acquire positive-valued imaginary parts. Furthermore, as L increases 
so the imaginary parts of the HCM parameters e^^^l'^ become increasingly large. Thus, we 
infer that the scattering interactions of increasing numbers of component phase inclusions 
become correlated as L increases, resulting in an increasingly dissipative HCM. The real parts 
of the HCM permittivity parameters e^j,|'^ do not vary significantly as the correlation length 
increases, as is illustrated in Figure 2. Notice that the relationship between the correlation 
length L and the second order SPFT parameters ej^^g'^ for rj = closely resembles the 
relationship between the inclusion size parameter rj and the zeroth order SPFT parameters 
^abo^ presented in Figure 1. 

The result of combining an inclusion size parameter rj > with a correlation length L > is 
presented in Figure 3. Therein, the components of the second order SPFT estimate e^j^ are 
plotted against both relative correlation length koL and the relative inclusion size 77 /l for 
the fixed volume fraction /„ = 0.3. The real and imaginary parts of the HCM permittivity 
parameters e^^|'^ are seen to increase with both 1] and L. Furthermore, the rate of increase of 
Im (e^^l'^) exceeds the corresponding rate of increase observed in Figure 1 when L = 0, and 
in Figure 2 when 77 = 0. Hence, by taking account of both nonzero depolarization volume 
and correlation length, more radiation is scattered out of the coherent field than would be 
the case if either nonzero depolarization volume or correlation length alone was taken into 
account. 



4.2 Dissipative component phases 

While it is instructive to consider nondissipative component phases in order to highlight 
the particular effects of depolarization volume and correlation length, in reality all mediums 
exhibit dissipation. Therefore, let us now consider the homogcnization of component phases 
characterized by the complex-valued permittivities = (2 + 0.05i) and ej, — (12 + OAi) Cq. 
We continue with the same values for the inclusion shape parameters Ux, Uy and Uz as given 
in (31). The volume fraction is fixed at /„ = 0.3 for all results presented in this section. 

In Figure 4, the components of the zeroth order SPFT estimate of HCM permittivity el'^j^ 
are plotted against relative inclusion size A:,,?]. Since the component phases are themselves 
dissipative, the imaginary parts of e^^g'^ are positive-valued at 77 = 0. As the inclusion size 
parameter rj increases, so the imaginary parts of e'^^j'^^ also increase. The growth in Im {^aW^) 
is attributable to the increasing amount of coherent scattering loss which develops as the 
depolarization regions increase in size. As was observed in the corresponding instance for 
nondissipative component phases, the real parts of e^^^ are relatively insensitive to 77. 

The second order SPFT estimates of the HCM permittivity parameters £^^2'^, calculated with 
7/ = 0, are graphed against relative correlation length A;„L in Figure 5. The observed increase 
in Im (c^^q'^) as L increases reflects the fact that coherent losses increase as the actions of 
increasing numbers of scattering centres become correlated. The real parts of the HCM 
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permittivity parameters vary little as the correlation length increases. Comparing Figures 4 
and 5, it is clear that the dependency upon L for the second-order SPFT with 77 = closely 
resembles the dependency upon r] for the zeroth order SPFT. 

Finally, let us consider the combined effect of nonzero depolarization volume and correlation 
length. In Figure 6, e^^|'^ are plotted against relative inclusion size k^^r| for the fixed relative 
correlation length k^L = 0.1. While Re (e^^|'^) are relatively insensitive to increasing 77, 
Im (ej^^l^) increase markedly as the inclusion size parameter increases. Furthermore, the rate 
of increase of Im (e^^|'^) in Figure 6 is greater than that seen in Figure 4 where the influence 
of 77 alone is considered, and in Figure 5 where the influence of L alone is considered. 

5 Concluding remarks 

In implementations of homogenization formalisms, such as the widely-used Maxwell Garnett 
and Bruggeman formalisms, the distributional statistics and sizes of the component phase 
particles are often inadequately taken into account. A notable exception is provided by the 
SPFT in which the distributional statistics of the component phases are described through a 
hierarchy of spatial correlation functions. In the present study the SPFT is extended through 
attributing a nonzero volume to the component phase inclusions. Thereby, we show that 
both depolarization volume and correlation length contribute to coherent scattering losses. 
Furthermore, the effects of depolarization volume and correlation length upon the imagi- 
nary parts of the HCM constitutive parameters are cumulative and reflect the anisotropy 
of the HCM. We remark that our numerical results which relate to the effect of depolar- 
ization volume alone, as presented in presented in Figures 1 and 4, are consistent with 
numerical results based on the extended Maxwell Garnett and extended Bruggeman homog- 
enization formalisms for isotropic dielectric HCMs [20, 32] . The importance of incorporating 
microstructural details within homogenization formalisms is thus further emphasized. 
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Figure 1: Real (left) and imaginary (right) parts of the zeroth order SPFT estimate of HCM 
permittivity e}^^ = eodiag(e^^Q, e^^^g, e^^g) plotted against volume fraction fa and relative 
inclusion size koT]. Component phase parameter values: = 2eo, eb = 12eo, = I, Uy = 3 
and Uz = 15. 
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Re£^t,2 




Figure 2: Real (left) and imaginary (right) parts of the second order SPFT estimate of HCM 
permittivity e^^j^ = eodiag(e^^25 ^afe2' ^ab2) plotted against volume fraction fa and relative 
correlation length k^L, with inclusion size parameter 77 = 0. Component phase permittivities 
ea and e^, and shape parameters Ux, Uy and f/^, as in Figure 1. 
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Figure 3: Real (left) and imaginary (right) parts of the second order SPFT estimate of 
HCM permittivity e'^j^ = eo diag (e^j,25 ^ab2' ^062) plotted relative correlation length k^L and 
relative inclusion size parameter rj/L, with volume fraction fa = 0.3. Component phase 
permittivities ea and eb, and shape parameters U^, Uy and Uz, as in Figure 1. 
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Figure 4: Real (top) and imaginary (bottom) parts of the zeroth order SPFT estimate of 
HCM permittivity e'^j^ = eodiag(e^^Q, e^j,Q, e^y^ plotted against relative inclusion size fco^ 
for volume fraction /„ = 0.3. Component phase parameter values: = (2 + 0.05i) eo, 
ej, = (12 + 0.4i) eo, = \, Uy = ?> and U.^ = 15. Key: e^^, e^^^ and elf^ are denoted by 
solid, broken dashed and dashed lines respectively. 
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Figure 5: Real (top) and imaginary (bottom) parts of the second order SPFT estimate of 
HCM permittivity e'^j^ = eodiag(e^^2) ^a62> ^062) plotted against relative correlation length 
koL, with inclusion size parameter 77 = and volume fraction fa = 0.3. Component phase 



permittivities and e?,, and shape parameters f/^, Uy and Uz, as in Figure 4. 
and e^^2 denoted by solid, broken dashed and dashed lines respectively. 
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Figure 6: Real (top) and imaginary (bottom) parts of the second order SPFT estimate of 
HCM permittivity e^j^ = eodiag(e^^2) ^avii ^062) plotted relative inclusion size parameter k^r]^ 
with relative correlation length k^L — 0.1 and volume fraction = 0.3. Component phase 
permittivities and e;,, and shape parameters [4, Uy and [4, as in Figure 4; and key as in 
Figure 5. 
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